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ABSTRACT 

We investigate the secular evolution of the orbital semi-major axis and ec- 
centricity due to mass transfer in eccentric binaries, allowing for both mass and 
angular momentum loss from the system. Adopting a delta function mass trans- 
fer rate at the periastron of the binary orbit, we find that, depending on the initial 
binary properties at the onset of mass transfer, the orbital semi-major axis and 
eccentricity can either increase or decrease at a rate linearly proportional to the 
magnitude of the mass transfer rate at periastron. The range of initial binary 
mass ratios and eccentricities that leads to increasing orbital semi-major axes 
and eccentricities broadens with increasing degrees of mass loss from the system 
and narrows with increasing orbital angular momentum loss from the binary. 
Comparison with tidal evolution timescales shows that the usual assumption of 
rapid circularization at the onset of mass transfer in eccentric binaries is not jus- 
tified, irrespective of the degree of systemic mass and angular momentum loss. 
This work extends our previous results for conservative mass transfer in eccentric 
binaries and can be incorporated into binary evolution and population synthesis 
codes to model non- conservative mass transfer in eccentric binaries. 

Subject headings: Celestial mechanics, Stars: Binaries: Close, Stars: Mass Loss 
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Introduction 



Many binary systems pass through at least one mass-transfer phase during the course 
of their evolution. These mass transfer episodes not only affect the internal evolution of 
the stellar components, but also impact the binary properties. In particular, changes in the 
mass ratio and transfer of linear and angular momenta between the stars cause changes in 
the orbital elements which affect the evolution of the binary and can cause a feedback on 
the mass transfer process. 

In binaries with eccentric orbits, the stars are closest to each other at the periastron of 
their relative orbit, so that any mass transfer is expected to take place first during periastron 
passage. The standard assumption in current binary evolution and population synthesis 
codes to deal with such mass transfer phases is that the orbit circularizes instantaneously at 
the onset of mass tran sfer. This assumption is in contrast with recent theoretical findings 



( jSepinsky et al. 



orbital eccentricities (IPetrova fc Orlov 



2007b) as well as observations of semi-detached binaries with non-zero 



1999 



Raguzova fc Popov 



2005|) 



Sepinsky et al. 



(l2007bl . hereafter Paper I) studied the orbital evolution due to mass 
transfer in eccentric binaries by deriving a set of perturbed equat ions of motion for the 



binary components, as outlined initially by 



Hadjidemetrioul (119691 ). The authors found 



that, under the assumption of conservation of total system mass and orbital angular 
momentum, the orbital semi-major axis and eccentricity can increase as well as decrease, 
depending on the initial orbital elements, binary component masses, and donor rotation 
rate. Furthermore, the orbital evolution timescales can be short enough to compete with or 
enhance any tidally driven orbital evolution. 

In this paper, we extend the analysis presented in Paper I to account for mass and 
angular momentum loss from the binary. In §[2] and [31 we briefly recall the relevant 
ingredients for the study of mass transfer in eccentric binaries derived in Paper I, and 
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update the formalism to account for systemic mass and angular momentum loss. In §U we 
present timescales for the evolution of the orbital semi-major axis and eccentricity due to 
mass transfer for different degrees of mass and angular momentum loss and we compare the 
timescales to the timescales of orbital evolution due to tidal dissipation. The final section 
is devoted to concluding remarks. 



2. Basic Assumptions 

We consider a close binary consisting of two stars with masses M\ and M 2 in an 
eccentric orbit with period P rb, semi-major axis a, and eccentricity e. The stars are 
assumed to rotate uniformly with angular velocities Q± and Q2 around an axis perpendicular 
to the orbital plane and in the same sense as the orbital motion. Since, in eccentric binaries, 
the magnitude of the orbital angular velocity f2 or b is a periodic function of time, the stellar 
rotation rates cannot be synchronized with the orbital motion at all orbital phases. 

At some time t, one of the stars is assumed to fill its Roche lobe initiating mass transfer 
to its companion through the inner Lagrangian point L\. We assume this point to lie on 
the line connecting the mass centers of the stars, even though non-synchronous rotation 
may cause i t to oscillate in the direction perpendicular to the orbital plane, as shown in the 



appendix of iMatese fc Whitmird (Il983l ). Since the donor's equatorial plane coincides with 
the plane of the orbit, the transferred mass can furthermore be assumed to remain in the 
orbital plane at all times. In what follows, we refer to the Roche lobe filling star as star 1 
and to the companion star as star 2. 

In Paper I, we assumed that all mass transfered by the donor was accreted by the 
companion, and that any orbital angular momentum transported by the transferred matter 
was immediately returned to the orbit (presumably through an accretion disk). Here, we 
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relax those assumptions and assume some fraction (3 of the transferred mass to be lost from 
the system: 

M T = (3M U (1) 

where Mt = Mi + M 2 is the total system mass. Consequently, the amount of matter 
accreted by the companion is 

M 2 = - 7 M! (2) 

where 7=1 — /?. The angular momentum carried away by the mass lost from the system is 
parameterized in terms of the specific angular momentum of the orbit as 

J^ = ^M T . (3) 



In the particular case where the matter lost from the system carries t 



Kolb et al 



re specific orbital 



200lh . For our purpose, 



angular momentum of the accretor, /1 = Mi/M 2 (e.g. 
we assume that no other sources of angular momentum loss besides mass loss are operating 
on the system. Any other sinks of orbital angular momentum such as tidal interactions, 
magnetic braking, and gravitational radiation can, at the lowest order of approximation, be 
added to Eq. (J3J) and Eqs. ffT8l) - fTT9l) derived in the next section to obtain the total rate of 
change of the orbital angular momentum and of the orbital elements. 

With these assumptions, the equations governing the motion of the two stars around 
their common center of mass can be written in the form of a perturbed two-body problem 

as 



d 2 f G (Mi + M. 



2, 



r + Si + Tjj+wi (4) 

where G is the Newtonian constant of gravitation, r is the position vector of the accretor 
with respect to the donor, x is a unit vector in the direction of r, y is a unit vector in 
the orbital plane perpendicular to r in the direction of the orbital motion, and z is a unit 
vector perpendicular to the orbital plane parallel to and in the same direction as Q rb- The 
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functions S, T, and W are the components of the perturbing force arising from the mass 
transfer between the binary components. In Paper I, we have shown that they are equal to 
c fl,x fl,x . M 2 ( ,=> , . 

- _t,,^ + — |r*|cos* — K|, (5) 

ry-i f2,y fl,y M 2 ( .=> , . , 

= m^~m7 + ^(^+1^11^1008. 

Mi / a Af 2l _ , 



2 



M 2 Mi' 1 J 

where Ai denotes the point on the donor's surface from which mass is lost (the L\ point), 

A 2 denotes the point on the accretor's surface at which mass is accreted, ta 1 and ?a 2 are 
the position vectors of A\ and A 2 with respect to the donor's and the accretor's center of 
mass, respectively, and the subscripts x, y, and z denote vector components in the x, y, 
and z directions. Moreover, f\ is the gravitational force exerted by particles in the mass 
transfer stream on the donor star, f 2 the gravitational force exerted by particles in the mass 
transfer stream on the accretor, VgMi the velocity of the matter ejected at Li with respect 
to the mass center of the donor star, vsm 2 the velocity of the accreted matter at A 2 with 
respect to the mass center of the accretor, and the angle between x and the vector from 
the center of mass of the accretor to A 2 . More details, as well as a diagram relating the 
above vectors, can be found in Paper I. 



3. Orbital Evolution Equations 

3.1. Secular Variation of the Orbital Elements 

The perturbing force with components S, T, and W in the equations governing 
the motion of the two stars around their common center of mass causes changes in 
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the orbital semi-major axis a anc 


eccen 


tricity e at rates 


Brouwer & Clemence 


1961; 


Danbv 


1962; 


Fitzpatrick 


1970) 



da 
~dl 



n(l-e 2 y/ 



Yj^[Se sin v + T(l + e cos z/)], 



(8) 



3 2U/2 



x <J 5 sin v + T 



2 cos v + e (1 + cos z/) 



(9) 



1 + e cos ^ 

where n = 27r/P or b is the mean motion, and v the true anomaly. The long-term secular 
evolution of the orbital elements is obtained by averaging these equations over one orbital 
period: 



da \ 
dt/ £ 

de \ 



p ^l 2 da 



-Porb J~P olh /2 dt 



dt. 



P orb /2 de 



-Porb J~P olb /2 dt 



dt. 



(10) 



As in Paper I, the integrals are most conveniently computed in terms of the true anomaly v. 



3.2. Orbital Angular Momentum Loss 

The perturbing functions S and T depend on the properties of the mass transfer 
stream. Calculation of the orbital semi-major axis and eccentricity evolution therefore, 
in principle, requires the calculation of the trajectories of the particles in the stream. In 
Paper I we bypassed such a calculation by assuming conservation of total system mass and 
orbital angular momentum. Here, we relax this assumption and generalize the formalism 
presented in Paper I by parameterizing systemic mass and angular momentum loss by 
means of Eqs. ([2]) and (J3]). 

As shown in Paper I, the rate of change of the orbital angular momentum is related to 
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the perturbing function T by 

-/orb _ Mi M 2 I Mr | (l-e 2 ) 1/2 y 

J orb Mi M 2 2M T na(l + ecosz/) ' 1 J 

If the only sink of orbital angular momentum is mass loss from the system, elimination of 
J or b from this equation and Eq. ([3]) yields 

^ 2 /Mr Mi M 2 na(l + e cosi/) ' 1 ' 



4. Orbital Evolution Timescales 

Since in binaries with eccentric orbits, mass transfer is expected to occur first at the 
periastron of the binary orbit, we approximate the mass transfer rate by a Dirac delta 
function as 

Mi = M 5(^), (14) 

where M < 0, |Mi is the instantaneous mass loss rate of star 1, and 5(u) the Dirac delta 
function. 

To calculate the rates of secular change of the orbital semi-major axis and eccentricity, 
we neglect the gravitational force exerted by the particles in the mass-transfer stream on 
the binary components and set 

fl,x = /2,* = 0, (15) 
fi,y = /2,y = 0. (16) 

Since the particles are confined to the orbital plane, the z-component of the perturbing 
forces are zero as well. A calculation of the actual gravitational effect of the particles in the 
mass transfer stream is beyond the scope of this paper and will be left for future studies. 
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With these assumptions, elimination of the perturbing function T between Eqs. 
and f[T3"j) . substituting Eqs. (j2J) and Q, and averaging over one orbital period yields 



iqV8M 2 ,y + v 5M 1; y 



11(1 



1/2 



|^orb,p| |rAi,p| 



7 g-l + (l- 7 ) (//+- 



7?l^orb,p||^4 2 | COS^p 



dv 



(17) 



where the subscript P indicates quantities evaluated at periastron, and q = M\jMi is the 
binary mass ratio. In the limiting case of conservative mass transfer (7 = 1), Eq. (117j) 
reduces to Eq. (35) in Paper I. 

Substituting Eqs. (|T5 i) -( fT7l) into Eqs. (jSJ) and (jBJ) yields expressions for the S and T 
components of the perturbing force acting on the binary due to the mass transfer between 
the binary components. With these expressions, Equations (TH|) (TTTj) yield the following rates 
of secular change of the orbital semi-major axis and eccentricity due to mass transfer at the 
periastron of eccentric binaries: 

aM 1 



da 
~dt 



ttM 1 (1 -e 2 )V2 



X 



a 



+ 7ge 



+ ( 7 g-l)(l- e 2 ) 



1 



COS $; 



q 



+ (i- 7) ^ +i )(i-^) i 



de 

Tit 



X 



iq- 



2tt 



a, 



cos $ P + 



\ r A u P\ 



+ 2( 79 -l)(l-e) 

+ 2(1-7)^+^(1-6)^ 



(19) 
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The rates of orbital evolution are thus linearly proportional to the magnitude Mq of the 
mass loss rate at periastron and dependent on the orbital semi-major axis a, the orbital 
eccentricity e, the donor mass M 1; and the binary mass ratio q. They also depend on 



the donor's rotational ang ular velocity Q 



Lagrangian point L\ (e.g. 



Sepinsky et al 



throu gh the position vector f f A 1 of the inner 



2007al ). Since I^a^pI oc a, the timescales explicitly 



depend on the orbital semi-major axis a only through the ratio jr^j/a of the accretor's 
equatorial radius to the orbital semi-major axis. In the limiting case of conservative mass 
transfer (7 = 1), Eqs. ffl8l) and ffl9|) reduce to Eqs. (39) and (40) in Paper I. 

In Figure [U the evolutionary timescales r a = a /a and r e = e/e for the semi-major 
axis a and orbital eccentricity e due to mass transfer in an eccentric binary are shown as a 
function of the binary mass ratio q, for different values of the orbital eccentricity and the 
parameters 7 and fi [see Eqs. 02) and (jBJ)]. While the actual timescales are given by the 
absolute values of r a and r e , we here allow the timescales to be negative as well as positive 
to distinguish between negative and positive rates of change of the orbital elements. 

For the calculation of the timescales, we assume the donor to rotate synchronously with 
the orbital angular velocity at periastron, and take the accretor to be a 1.44 M & neutron 
star. The donor mass is then fixed by the binary mass ratio. Since the radius of the neutron 
star is much smaller than the semi-major axis of the orbit, the terms in Eqs. (1181) and (fT9l) 
containing the ratio \rA 2 \/ a are negligible compared to the other terms. We therefore set 
\ta 2 \ = 0, so that the timescales r a and r e are independent of a. For the mass transfer rate 
we adopt a constant Mq = — 10~ 9 M Q yr _1 , but note that the linear dependence of the rates 
of change of the orbital semi-major axis and eccentricity on M allows for an easy rescaling 
of the timescales to different mass transfer rates. 

The overall shape of the curves shown in Figure [1] is similar to that of the curves shown 
in Figure 2 of Paper I in the case of conservative mass transfer. For a given fraction 7 of 
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Fig. 1. — Orbital evolution timescales for the semi-major axis a (top) and orbital eccentricity 
e (bottom) for a delta function mass loss rate M\ = M 5(u) with M = — lO _9 M yr _1 
for a range of eccentricities e = 0.1, 0.3, 0.5, 0.7, 0.9, as shown in the bottom middle 
panel of the left plot. The timescales are calculated as a function of the binary mass ratio 
q = M1/M2 under the assumptions that the donor star rotates synchronously with the 
orbital angular velocity at periastron and that the accretor is a neutron star of mass M 2 = 
1.44 Mq. The different panels show timescales for different values of 7 and /i. The \i — q 
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mass loss from the system, increasing /i implies more angular momentum loss from the 
binary, and thus faster orbital shrinkage or slower orbital expansion. For a given degree \i 
of specific orbital angular momentum loss, on the other hand, decreasing 7 implies more 
mass loss from the binary, causing faster orbital expansion and slower orbital contraction. 
The timescales for the evolution of the orbital eccentricity show a similar dependence on 
the parameters 7 and \i. 

The evolutionary timescales of the orbital semi-major axis a and eccentricity e can be 
negative as well as positive, depending on the initial binary properties. The transition from 
shrinking to growing orbital elements is illustrated more clearly by the contour plots shown 
in Figure [21 In these plots, the gray shades represent different timescales of orbital evolution 
and the thick black line marks the transition from shrinking [right of the thick black line) to 
growing (left of the thick black line) orbital elements. For comparison, the dashed black line 
marks the transition in the case of fully conservative mass transfer presented in Paper 3 
As 7 decreases, more mass is lost from the system and the transition from negative to 
positive rates of change of the semi- major axis and eccentricity moves to larger mass 
ratios. Conversely, as fi increases, more angular momentum is lost from the system and the 
transition from negative to positive rates of change of the semi-major axis and eccentricity 
moves to smaller mass ratios. In the particular case where matter leaving the system carries 
away the specific orbital angular momentum of the accretor (fi = q), the critical mass ratio 
separating increasing from decreasing orbital semi-major axes and eccentricities increases 



x The sharp bend of the dashed line toward smaller q values for e < 0.95 was not observed 
in Paper I due to the lower resolution of the orbital eccentricity grid considered in that paper. 
Equations (41) and (42) of Paper I still fit the line to better than 1.5% (10%) for e < 0.8 
(e < 0.95) in the case of the orbital semi-major axis and for e < 0.85 (e < 0.95) in the case 
of the orbital eccentricity. 
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Fig. 2. — Contour plots of the orbital evolution timescales for the semi-major axis a (top) 
and orbital eccentricity e (bottom) in the (q, e)-plane for the same set of assumptions as 
adopted in Fig. [TJ From the darkest to the lightest gray, the different gray shades represent 
timescales from to 1 Gyr, 1 to 5Gyr, 5 to 10 Gyr, 10 to 15 Gyr, and more than 15Gyr, 
respectively. The thick black line in each panel separates the regions of the (q, e) space where 
a > or e > (left of the thick black line) from the regions of the (q, e) space where a < 
or e < (right of the thick black line). The timescales in the bottom left panel (7 = and 
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by about 30% when going from fully conservative to fully non-conservative mass transfer. 

To assess the role of mass transfer in the overall evolution of the binary, we compare 
the timescales shown in Figures [1] and [2] with the orbital evolution timescales due to tidal 



di ssipation in t 



in 



r e Roc he-lobe filling star. The tidal evolution timescales are determined as 



Hurley et al. 



(120021 ) assuming the donor is a zero-age main-sequence star. The variations 
of the tidal evolution timescales as a function of the orbital eccentricity and binary mass 
ratio are shown in Figure 4 of Paper I. The discontinuity of the timescales at q ~ 0.87 
corresponds to the transition from donor stars in which convective damping is the dominant 
tidal dissipation mechanism (M 2 < 1.25 M Q ) to donor stars in which radiative damping is 
the dominant tidal dissipation mechanism (M 2 > 1.25 M Q ). 

The timescales of orbital evolution due to the combined effects of mass transfer and 
tides are shown as contour plots in Figure [3j As before, the thick black line indicates 
the transition from shrinking (right of the thick black line) to growing (left of the thick 
black line) orbital semi-major axis and eccentricity. For comparison, the transition line for 
conservative mass transfer is shown by means of the dashed black line. The vertical white 
line at q ~ 0.87 separates donor stars in which tidal energy is dissipated by convective 
damping from donor stars in which tidal energy is dissipated by radiative damping. As 
noted in Paper I, in the case of conservative mass transfer, there are large regions in the 
(q, e) parameter space where the combined effects of mass transfer and tidal interactions 
do not lead to rapid circularization of the orbit after the onset of mass transfer and where 
eccentricity pumping occurs instead of eccentricity damping. When mass loss from the 
system is taken into account, the parameter space for eccentricity pumping becomes even 
larger, though the timescales of orbital evolution in the newly accessible e > regions are 
long (> 15Gyr). Similar behavior is observed for the evolution of the orbital semi-major 
axis. 
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Fig. 3. — Contour plots of the orbital evolution timescales for the semi-major axis a (top) 
and orbital eccentricity e (bottom) in the (q, e)-plane due to the combined effects of mass 
transfer and tidal dissipation in the donor star of a semi-detached binary. The mass transfer 
rate is assumed to be a delta function of amplitude Mq = — 10 -9 M Q yr -1 at the periastron of 
the binary orbit, the donor star is assumed to rotate synchronously with the orbital angular 
velocity at periastron, and the accretor is assumed to be a 1.44 M neutron star. The tidal 
contribution to the orbital evolution timescales is determined assuming the donor star has 
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We note that the results for zero-age main- sequence stars presented here are intended 
for illustration and comparison with Paper I. As tides will be much more efficiently damped 
in giant stars with deep convective envelopes, the interplay between mass transfer and tidal 
interactions will likely change significantly with the evolutionary state of the donor star. 



5. Concluding Remarks 



We extended the formalism to study the orbital evolution due to mass transfer 
in eccentric binaries derived in Paper I to account for the effects of mass and angular 
momentum loss from the system. Adopting a delta function mass transfer rate at the 
periastron of the binary orbit, we find that the usually adopted assumption of rapid 
orbital circularization during the early stages of mass transfer remains unjustified when 
systemic mass and angular momentum loss are taken into account. Our results thus 
present a possible explanation for the observation of non-zero orbital eccentricities in 
mass-transferring binaries. 

The formalism presented in this paper and in Paper I can be incorporated into binary 
evolution and population synthesis code to provide a model for eccentric mass-transferring 
binaries which are currently, by construction, absent in any population synthesis studies 
of interacting binaries and their descendants. In future work, we intend to incorporate 
this orbital evolution into such code, and determine the prevalence and the long-term 
evolutionary effects of mass transfer in eccentric binary systems. Possible applications 
include the modeling of systems such as CirX-1 in which a neutron star is thought to 
accrete matter from a Roche-lobe filling or nearly Roche-lobe filling companion during 



each periastron passage and in which near-IR anc 



of an accretion-driven mass outflow f Clark et al. 



X-ray spectroscopy support the prese nce 



2003 



Iaria et al. 



2008 



Tennant et al 



19861 ). Another example is the ultracompact X-ray binary 4U 1820-30 which is thought 



to be a member of a hierarchical triple (jChou fc Grindlay 



2001 



Zdziarski et al. 



20071 ) 



Furthermore, dynamical interactions between single and binary stars in dense stellar clusters 
can indu ce eccentricities in circ ular binaries and enhance eccentricities in already eccentric 



binaries (IHeggie fc Rasio 



19961 ). This induced eccentricity can directly lead to Roche Lobe 



overflow at perias tron, as has been suggested for the flaring X-ray binaries in NGC 4697 



( iMaccarone 



20051 ). In future work, we intend to model the mass transfer rate at periastron 
more realistically by taking into account the atmospheric properties of the donor star and 
considering the feedback of the orbital and radial evolution of the star on the mass transfer 



rate. 
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